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Within the context of the Abelian Projection of QCD monopole-like quantum excitations of 
gauge fields are studied. We start with certain classical solutions, of the SU(2) Yang-Mills field 
equations, which are not monopole-like and whose energy density diverges as r — > do. These 
divergent classical solutions are then quantized using a modified version of Heisenberg's quantization 
technique for strongly interacting, nonlinear fields. The modified Heisenberg quantization technique 
leads to a system of equations with mixed quantum and classical degrees of freedom. By applying 
a Feynman path integration over the quantum degrees of freedom the quantum-averaged solution 
' gives a nondivergent, monopole-like configuration after Abelian Projection. 
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Monopoles have been studied within both Abelian and non- Abelian @ gauge theories. Usually one starts with a 
classical monopole configuration (i.e. having a magnetic field which becomes Coulombic as r — > oo) and then consider 
the quantum corrections to the system. One difficulty with this approach is that the coupling strength of monopole 
configurations is large due to the Dirac quantization condition Q which requires that there be an inverse relationship 
between electric and magnetic couplings. A perturbatively small electric coupling requires a non-perturbatively large 
t-H i magnetic coupling. The perturbative quantization techniques do not work well with monopoles for the same reason 
that they have trouble with QCD in the low energy regime : the couplings are non-perturbatively large. 

Numerical simulations in lattice QCD show that the configurations of gauge fields with monopoles dominate in 
the path integral, i.e. monopole configurations are important contributors to the QCD string tension. Thus the 
expectation value of any physical quantity (T) in non-Abelian gauge theory can be accurately calculated (see for 
review Ref. ||) in terms of the monopole currents extracted from the projected Abelian fields. Magnetic monopoles 
necessarily emerge as relevant degrees of freedom in the Abelian Projections of lattice gluodynamics. In this paper 
we present a specific picture of how such magnetic monopole configurations can emerge and dominate the QCD path 
integral. We start with spherically symmetric classical solutions to the SU(2) field equations. These solutions have 
a divergent energy density. However in constrast to the Wu-Yang monopole solutions, which are divergent due to 
fields that blow up at r — 0, our solutions are divergent due to fields which blow up at r = oo. These divergent 
^ , solutions are parameterized by two ansatz functions. One of the ansatz functions smoothly and monotonically diverges 
at oo while the other functions strongly oscillates. Applying a modified variant of a non-perturbative quantization 
method originally used by Heisenberg Q 0], we find that the quantized version of our divergent classical solution 
becomes well behaved: the divergent ansatz function now goes to at oo, and the rapid ocsillations of the other ansatz 
function are "smoothed" out. Rather than quantizing an SU(2) field configuration which is monopole-like already 
at the classical level, we quantized (using Heisenberg's method) a non-monopole configuration, which has fields and 
an energy density which diverge at spacial infinity (r — > oo). After quantization this classically divergent solution 
becomes physically well behaved, and its asymptotic magnetic field becomes monopole-like. This may indicate that 
monopoles are inherently quantum objects: rather than quantizing a configuration which is monopole-like at the 
classical level, real monopoles might arise as a consequence of quantizing non-monopole classical solutions. 



II. CLASSICAL SU(2) SOLUTION 

In this section we will briefly review the classical, spherically symmetric, SU(2) Yang-Mills theory solutions to 
which we will apply the modified Heisenberg quantization method. We begin with the following ansatz for pure SU(2) 
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Yang-Mills theory 

K = ^ 9 {r), (1) 

K = e aij ^[l-f(r)], (2) 

here a = 1,2,3 is a color index; the Latin indices k,l — 1, 2, 3 are space indices. This is the Wu-Yang ansatz |J. 
Under this ansatz the Yang - Mills equations (D^F a M „ = 0) become 

r 2 f" = f-f-fg 2 , (3) 

r 2 g" = 2gf (4) 

where the primes indicate differentiation with respect to r. In the asymptotic limit r — > oo the general solution to 
Eqs. (^), ([|) approaches the form 

f{x) « Asin^+^o), (5) 

A 2 = a(a - 1) (7) 

where x = r/rg is a dimensionless radius and fo, 0o> an d A are constants. The plot of these functions are presented 
on the Fig.g |. 




FIG. 2: Function g(x) 
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The numerical study of Eqs. (||), (^) showed that the approximate form of the solution in Eqs. (J5j) - (Q) was a good 
representation even for r close to the origin. These solutions can be compared with the certain confining potentials 
used in phenomenological studies of quarkonia bound states fij . It is possible to find analytic solutions with increasing 
gauge potentials, similar to the above numerical solution for the ansatz g(x), if one couples the Yang-Mills fields to a 
scalar field J^]. However, it is the increasing of g(r) that leads to the field energy and action of this solution diverging 
as r — > oo. 

The "magnetic" (Hf = tijkFjj.) an< ^ "electric" (Ef = FqA fields associated with this solution can be found from 
the non-Abelian gauge potentials, A^, and have the following behavior 

H r* f —^-i H«ocf, H$Kf', (8) 
E^ !X 7 t^l > E^oc^-, ESx^l. (9) 



III. EXTRACTION OF THE OFF-DIAGONAL COMPONENTS OF GAUGE POTENTIAL 

Our SU(2) gauge potential (P (||) can be written as 

K = [l-/(r)] {ship; -cos *>;()}, (10) 
A q v = [1 — f(r)] sinO {cos 9 cos ip; cos 6 simp; — sind} , (11) 



r 



The gauge transformation 



with the matrix 



leads to the following result 



A'l = S- 1 AlS-i(d ll S- 1 )S (13) 

(14) 



cos § —e* sm ■ 
e iv sin | cos | 



X = ( „ ( l " () ■ U- r >> 



^ = (costf - 1) r _ U J - f(r) sine (V n ) , (16) 
ff( r ) A 



4 - r [o -l) ■ <") 

In this form we see that if f(r) — > as r — > oo then only the first term in Eq. ( |l6| ) is relevant. This first term is in 
the U(l) subgroup of the SU(2) gauge field, and gives rise to a magnetic monopole. 

Under the assumption (originally found in ref. JjJ) that the U(l) diagonal components of the non-Abelian gauge field 
are the most important for confinement, one tries to make a gauge choice which minimizes or excludes the off-diagonal 
components of the gauge potential (i.e. the Cartan subalgebra). For example, in Maximal Abelian Projection [flo|| , 
JTl| the following gauge condition is applied 

D, (A H ) A°Jf = d»A°Jf + ig [<, A°/f] = (18) 

where A^ = A^ + A°ff . A^ is the Abelian, diagonal pomponent of A^ while A^ff are the off-diagonal components 
of A p. This gauge corresponds to the minimization of the following functional 

R 



J {K] 2 + K] 2 } di - (is) 



where A°Jf — A 1 ^ 2 . Lattice calculations (for review, see for example, |Q) indicate that monopole- like gauge field 
configurations are dominate in the path integral. Here we show that one can reach this same conclusion by using a 
variant of Heisenberg's quantization method on non-monopole, classical solution. By quantizing the non-Abelian we 
find that f(r) — > as r — > oo. Thus the quantized non-Abelian fields asymptotically appear as monopoles. 
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IV. HEISENBERG'S NON-PERTURBATIVE QUANTIZATION 

Although the confining behavior of these classical solutions is interesting due to their similarity with certain phe- 
nomenological, confining potentials, the infinite field energy makes their physical importance/meaning uncertain. One 
possible resolution to the divergent field energy is if quantum effects removed the bad long distance behavior. The 
difficulty is that strongly interacting, nonlinear theories are notoriously hard to quantize. In order to take into account 
the quantum effects on these solutions we will employ a variation of the method used by Heisenberg 0] in attempts 
to quantize the nonlinear Dirac equation. It can be shown ]T^ | that the Gor'kov's equations (and consequently the 
Ginzburg - Landau equations) in superconductivity theory is the direct application of Heisenberg's idea for quantiza- 
tion of Cooper pairs in the superconductor. We will outline the key points of the method by using the nonlinear Dirac 
equation as an illustrative example. The nonlinear spinor field equation considered by Heisenberg had the following 
form : 

7 ^(x) - / 2 3#(^)] = (20) 

where 7 M are Dirac matrices; ip(x), ?/> are the spinor field and its adjoint respectively; ^s[ip(?p?p)] is the general nonlinear 
spinor self interaction term which involved three spinor fields and various combinations of j^'s and/or 7 5 's. The 
constant I has units of length, and sets the scale for the strength of the interaction. Next one defines r functions as 

T( Xl x 2 ...\ yi y 2 ...) = (o\T$(xi)$(x 2 )...ip*(yi)4>*(y2)...m (21) 

where T is the time ordering operator; $} is a state for the system described by Eq. (po|). Applying Eq. ( |2C| ) to j2l| ) 
we obtain the following infinite system of equations for various r's 

Q 

l~ tf r )-7yir T(xi...x n \y 1 ...y n ) = %{T(x 1 ...x n x r \y 1 ...y n y r )} + 

S(x r - y 1 )T(x 1 ...x r _ 1 x r+1 ...x n \y 2 --y r -iyr+i---yn) + 

S(x r - y 2 )T(xi...x r -ix r+ i...x n \yiy 2 ...y r ^iy r+ i...y n ) + ... (22) 

Eq. ( p2| ) represents one of an infinite set of coupled equations which relate various order (given by the index n) of the 
t functions to one another. To make some head way toward solving this infinite set of equations Heisenberg employed 
the Tamm-Dankoff method whereby he only considered r functions up to a certain order. This effectively turned the 
infinite set of coupled equations into a finite set of coupled equations. 

For the SU(2) Yang-Mills theory this idea leads to the following Yang-Mills equations for the quantized SU(2) gauge 
field 

D ^pa^ = q (23) 

here F a ^ v is the field operator of the SU(2) gauge field. 

One can show that Heisenberg's method is equivalent to the Dyson-Schwinger system of equations for small coupling 
constants. One can also make a comparison between the Heisenberg method and the standard Feynman diagram 
technique. With the Feynman diagram method quantum corrections to physical quantities are given in terms of 
an infinite number of higher order, loop diagrams. In practice one takes only a finite number of diagrams into 
account when calculating the quantum correction to some physical quantity. This standard diagrammatic method 
requires a small expansion parameter (the coupling constant), and thus does not work for strongly coupled theories. 
The Heisenberg method was intended for strongly coupled, nonlinear theories, and we will apply a variation of this 
method to the classical solution discussed in the last section. 

We will consider a variation of Heisenberg's quantization method for the present non-Abelian equations by making 
the following assumptions [ O : 

1. The physical degrees of freedom relevant for studying the above classical solution (which after quantization 
will be spherically symmetric excitations in QCD vacuum) are given entirely by the two ansatz functions /, g 
appearing in Eqs. (|3|), (^J). No other degrees of freedom will arise through the quantization process. 

2. From Eqs. (§), (§) we see that one function f(r) is a smoothly varying function for large r, while another 
function, g(r), is strongly oscillating. Thus we take g(r) to be an almost classical degree of freedom while f(r) 
is treated as a fully quantum mechanical degree of freedom. Naively one might expect that only the behavior 
of second function would change while the first function stayed the same. However since both functions are 
interrelated through the nonlinear nature of the field equations we find that both functions are modified. 
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To begin we replace the ansatz functions by operators f(x),g(x). 



r 2 f" = f-f-fg 2 , 
r g = 2gf 



(24) 
(25) 



These equations can be seen as an approximation of the quantized SU(2) Yang-Mills field equations (|22|), Taking into 
account assumption (2) we let g —* g so that it becomes just a classical function again, and replace f 2 in Eq. (25) by 
its expectation value to arrive at 



r 2 f" = f-f-fg 2 , 

r 2 9" = 2g(f) 



(26) 
(27) 



where the expectation value (f 2 ) is taken with respect to some quantum state \q) - (f 2 ) — (q\f 2 \q). We can average 
Eq. <M) to get 



2 (/)" 

r 2 g" 



(f) (/} 

2.9</ 2 } 



(f)g 2 



(28) 
(29) 



Eqs. (28) ( p9| ) are almost a closed system for determining (/) except for the (f 2 ) and (/ 3 ) terms. One can obtain 
differential equations for these expectation values by applying r 2 d/dr to f 2 or f 3 and using Eqs. (^6|) - (|27j). However 
the differential equations for (f 2 ) or (/ 3 ) would involve yet higher powers of / thus generating an infinite number of 
coupled differential equations for the various (/"). In the next section we will use a path integral inspired method 
p5| to cut this progression off at some finite number of differential equations. 



PATH INTEGRATION OVER CLASSICAL SOLUTIONS 



Within the path integral method the expectation value of some field $ is given by 



(30) 



The classical solutions (if they exist), <& c ;, give the dominate contribution to the path integral. For a single classical 
solution one can approximate the path integral as 



(31) 



where A is a normalization constant. Consequently the expectation of the field can be approximated by 

$ e iS[#].D$ w $ c/ . (32) 



(33) 



We are interested in the case where $ is the gauge potential A", in which case our approximation becomes 



4" 



iS 



00 



D 



■to 



(4)*, 



(A")^ are the classical solutions of the Yang - Mills equations labelled by a parameter (f> . In the present case the 
classical solution with the asymptotic form (||) (JsJ) has an infinite energy and action. When one considers the Euclidean 
version of the path integral above the exponential factor in ( |33| ) becomes exp[— SKA^)^] which for an infinite action 
would naively imply that this classical configuration would not contribute to the path integral at all. However, there 
are examples where infinite action, classical solutions have been hypothesized to play a significant role in the path 
integral. The most well known example of this is the meron solution which has an infinite action. Analytically 
the singularities of the meron solutions can be dealt with by replacing the regions that contain the singularities by 
instanton solutions. Since instanton solutions have finite action this patched together solution of meron plus instanton 
has finite action. The drawback is that the Yang-Mills field equations are not satisfied at the boundary where the 
meron and instanton solutions are "sewn" together. In addition recent lattice studies |f7j jl8| have indicated that 
merons (or the patched meron/instanton) do play a role in the path integral. Here we will treat this divergence in 
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the action in an approximate way through a redefinition of the path integral integration measure. We now define the 
path integral by the following approximate way 



(A a Jr)) = / Ale- S ^DA 



fe-^ L i( A ^ r) U dV d(A-(r) 



(34) 



where dV is the 3D infinitesimal volume. Here our basic assumption is that in the first rough approximation (A"(r)) 
can be approximately estimated by integration over classical singular solutions. Every such solution is labelled by <fio 
and consequently we should integrate over phase <f>Q. It can be shown that the Lagrangian for the system is 



L w a 2 (a - I) \ a - 1 + 2a sin 



„2a 



(35) 



The next our assumption is that we neglect the oscillating term sin 2 [(r/ro) a 
cancelled in numerator and denominator. That leads to 



i>o] and factor e 



-if a 2 (a-l) 2 r 2a dV 



Jdcbo (I«(r)) 



4>0 



2tt 



(36) 



here the measure d (^A^r^j 
various powers of f(r) 



dcftQ. Using this and the fact that A^ (r) ex /(r) we calculate the expectation value of 



(/) 
</ 2 > 
</ 3 > 
(/ 4 ) 



1 

2^ 
I 

2^ 



271 



cia(po 



d u 90 



A 
2^ 

£_ 

2n 



sin(x" + <f)o)d(f)o = 0, 



2 71 



sin (x a + (f>o)d(j) 



2 ' 



271 





271 



27 ' /cW ' n 



271 



A 3 

— / sin 3 (x a + o )d0 o = 0, 

27T 

.4 4 



o 

271 



i J / c 4 z #o = ^ I sin 4 (x« + 0oWo = lA 4 



(37) 
(38) 
(39) 
(40) 



where f c i is the function from the Eq.(^|). The path integral inspired Eqns. (|37]) - (|40|) are the heart of the cutoff 



procedure that we wish to apply to Eqns. 



28D 



On substituting Eqs. (|37|), (|3j) into Eqs. (g9|) we find 



that Eq.(|2g) is satisfied identically and Eq.(|29j) takes the form 

r 2 g" = A 2 g = a(a - l)g, 

which has the solutions 



a > 1 



fJ 
fJ 



9or , 
9o 



(41) 

(42) 
(43) 



where go is some constant. The first solution is simply the classically averaged singular solution (g) which still has 
the bad asymptotic divergence of the fields and energy density. The second solution, (43), is more physically relevant 
since it leads to asymptotic fields which are well behaved. 

The solution of Eq. (^) implies the following important result: the quantum fluctuations of the strongly 
oscillating, nonlinear fields leads to an improvement of the bad asymptotic behavior of these nonlinear 
fields. After quantization the monotonically growing and strongly oscillating components of the gauge potential 
became asymptotically well behaved. We interpret this spherically symmetric quantized field distribution as a model 
of the monopole excitations in QCD vacuum. 
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As r — > oo we find the following SU(2) color fields 

(H^ g ) cx (/') « 0, (45) 



H r) « — -J with 0= ga(a -1)-1, (44) 



(^r> cx (^^)«-^r, (46) 

(*V,9> « ^ = 0. (47) 

We can see that as r — > oo <C |(.ff")|. In particular at infinity we find only a monopole "magnetic" field 

H" « Q/r 2 with a "magnetic" charge Q. This result can be summarized as: £/ie approximate quantization 
of the SU(2) gauge field (by averaging over the classical singular solutions) gives a monopole-like 
configuration from an initial classical configuration which was not monopole-like. We will call this a 
"quantum monopole" to distinguish it from field configurations which are monopole-like already in the classical 
theory. 



VI. ENERGY DENSITY 



The divergence of the fields of the classical solution given by Eqs. (^|) - leads to a diverging energy density for 
the solution, and thus an infinite total field energy. The energy density e of the quantized solution is 



w2 , , ua ,2„(r9'-9Y , 2</V , 2</' 2 ) , ((/ 2 -I) 2 ) 



e «(^ + (fr2r«(-^) + + + ~ 4 ( 48 ) 



The first two terms on the right hand side of Eq. (ff8|), which involve the "classical" ansatz function g(r), go to zero 
faster than 1/r 4 as r — -> oo due to the form of g(r) in Eq. (^3|). Thus the leading behavior of e is given by the last 
two terms in Eq. (jig ) which have only the "quantum" ansatz function f(r). To calculate (J' 2 } let us consider 

rf^ (/'(r 2 )/'(r0) - </'(r 2 )/ 3 (r 1 )) - (/'(r a )/(n)) (l + g 2 (n)) (49) 

In the limit r' — > r we have 

(/'(r)/"(r)> = ^ (/'>))' = 1 (/ 4 (r))' - \ (/ 2 (r))' (l + , 9 2 (r)) . (50) 

From Eqs. ( |38| ) and ( ^0| ) we see that (f 2 )' = which implies (f 2 ) = const. Thus the third term in fl48| ) gives the 
leading asymptotic behavior as r — > oo to be 

const 

£ ~— ( 51 ) 

and the total energy of this "quantum monopole" (excitation) is infinite. This fact indicates that our approximation 
( |30| ) is good only for (f n (r)) calculations but not for the derivative (/' (r)}. 



VII. CONCLUSIONS 



Starting from an infinite energy, classical solution to the SU(2) Yang-Mills field equations we found that the bad 
asymptotic behavior of this solution was favorably modified by a variation of the quantization method proposed by 
Heisenberg to deal with strongly coupled, nonlinear field theories. In addition, although the original classical solution 
was not monopole-like, it was found that the quantized solution was monopole-like. One possible application of 
this is to the dual-superconductor picture of the QCD vacuum. In this picture one models the QCD vacuum as a 
stochastic gas of appearing/disappearing monopoles and antimonopoles as in Fig.||. These monopole/ antimonopole 
fluctuations can form pairsf analogous to Cooper pairs in real superconductors) which can Bose condense leading 
to a dual Meissner effect |i~9| [^(J expelling color electric flux from the QCD vacuum, expect in narrow flux tubes 
which connect and confine the quarks. Lattice calculations confirm such a model pl] j : monopoles appear to play a 
major role in the QCD lattice gauge path integral. Based on the results of the present paper it may be that : (a) the 
monopoles which are considered in the dual superconductor QCD vacuum picture and (b) the monopoles in lattice 
simulations whose location is determined by the Abelian Projection procedure should be the "quantum" monopoles 
discussed here. 
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